New superconducting and semiconducting Fe-B compounds predicted with an ab initio 

evolutionary search 
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New candidate ground states at 1:4, 1:2, and 1:1 compositions are identified in the well known 
Fe-B system via a combination of ab initio high-throughput and evolutionary searches. We show 
that the proposed oP12-FeB2 stabilizes by a break up of 2D boron layers into ID chains while oPlO- 
FeB4 stabilizes by a distortion of a 3D boron network. The uniqueness of these configurations gives 
rise to a set of remarkable properties: oP12-FeB2 is expected to be the first semiconducting metal 
diboride and oP10-FeB4 is shown to have the potential for phonon-mediated superconductivity with 
a T c of 15-20 K. 
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A range of advanced compound prediction methods 
has been developed recently to accelerate the experimen- 
tal search for materials displaying novel physics or tech- 
nologically relevant features [l|-|j]. Unconstrained struc- 
tural optimization with evolutionary algorithms (EAs) 
has shown the ability to predict complex configurations 
given only the composition leading to identification of 
exotic high-pressure phases Q . High-throughput screen- 
ing with data mining techniques has proven effective in 
revealing compositions favorable to form in large sets of 
multi-component systems In this study we demon- 
strate that new ambient-pressure materials with appeal- 
ing properties could be found in such a well-known and 
accessible binary system as Fe-B. 

The experimental research on Fe-B compounds has 
been driven primarily by their potential to serve as a 
hardening agent in steels [H or as hard protective coat- 
ings 0, 0. According to the latest experimental phase 
diagram |8j, FeB and Fe2B are the only reproducible low 
temperature phases that have been shown to crystallize 
in the oP8 (or the related oS8 Q) and tI12 configurations, 
respectively. Less is known about the boron-rich ordered 
phases with only a few reports available: observation of a 
metastable FeB 4 g intercalation compound [Io[ and pos- 
sible synthesis of amorphous [ll[ and the AlB 2 -type [l2[ 
iron diborides. Previous modelling work on Fe-B com- 
pounds has given insights into their binding, magnetic, 
and structural properties but has not systemati- 

cally explored the possibility of obtaining new stable iron 
borides. 

Our reexamination of the Fe-B system within density 
functional theory (DFT) begins with a high-throughput 
scan of known configurations listed in the Inorganic Crys- 
tal Structure Database ICSD [lij]. We show that never 
observed oP6-FeB 2 (IiP6-FeB 2 ) and tI16-FeB phases are 
marginally stable relative to the known compounds. The 
proposed Fe-B ground states are then refined with an ab 
initio evolutionary search that suggests oP10-FeB 4 and 
oP12-FeB 2 to be ground states at 1:4 and 1:2 composi- 



tions. The prediction of the brand new stable structure 
types is surprising as transition metal (TM) borides tend 
to crystallize in configurations correlating well between 
the 3<i, Ad, and 5d series [2(| . We link the stabilization of 
the Fe-B phases to the structural changes in the B net- 
works that lead to radically new properties. At 1:2 metal- 
boron composition, famous for the outstanding MgB 2 
superconductor 2l( and the hardest metal-based ReB 2 
material [22j, oP12-FcB 2 stands out as the first semicon- 



ducting metal diboride made out of B chains rather than 
B layers. At 1:4 composition, the nonmagnetic oPlO- 
FeB 4 is examined using electron-phonon (e-ph) calcula- 
tions and predicted to be, subject to spin fluctuation ef- 
fects [231 ] . a superconductor with an unexpectedly high 
T c of 15-20 K. The critical temperature falls between the 
typical 10 K T c of TM borides [H| and the 39 K T c of 
MgB 2 [2l|. If synthesized, 0PIO-FCB4 could extend the 
family of recently discovered iron-based LaFeAsOi_ K F x 
and FeSe superconducting materials [25[ but have the 
conventional phonon-mediated coupling mechanism. 

We carry out the high-throughput scan by calculating 
formation enthalpies at T ~ K and P = GPa with 
VASP [26[ for over 40 commonly observed M-B and M- 
C ICSD structure types [27| in the whole composition 
range. The B-rich end of the phase diagram is further 
explored with the Module for Ab Initio Structure Evo- 
lution (MA1SE) linked with VASP which enables an 
EA search for the lowest enthalpy ordered phases. The 
unconstrained structural optimization is carried out for 
most likely to occur 1:6, f:4, f:3, and 1:2 compositions 
starting from random unit cells of up to 15 atoms (for 
further details see supplementary material [13]). 

Finding ground states also depends on the accuracy 
of the simulation method and the inclusion of important 
Gibbs energy contributions [l|. We use the projector 
augmented waves method [29( and allow spin polariza- 
tion unless stated otherwise; the chosen ener gy c utoff of 
500 eV and dense Monkhorst-Pack fc-meshes [30j ensure 
numerical convergence of formation energy differences 
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FIG. 1: (Color online) Stability of Fe-B alloys calculated with 
GGA-PBE: a) formation enthalpy; b-d) Gibbs energy with 
thermodynamic corrections due to the vibrational entropy for 
selected candidate phases w.r.t. the a-B-H-oP8-FeB tie-line. 
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FIG. 2: (Color online) Competing B-rich Fe-B phases; cell 
parameters are given in the supplementary material [27j|. 



to typically 1-2 mcV/atom. Wc employ the Perdew- 
Burkc-Ernzcrhof (PBE) exchange-correlation (xc) func- 
tional [3ll ] within the generalized gradient approxima- 
tion (GGA) that provides a realistic description of the 
Fe grou nd state [321 ] . Tests in the supplementary mate- 
rial j27| demonstrate independence of our key finding, the 
stability of new phases at the 1:2 and 1:4 Fe-B compo- 
sitions with respect to known compounds, on the choice 
of the xc functional. The ground state of B is modelled 
as a-B which has been recently shown to be only 3-4 
meV/atom above the more complex /3-B in the 0-300 K 
temperature range [Hj]. We include phonon corrections 
to G(T) using a finite displacement method as imple- 
mented in PHON [3^ . The e-ph calculations are carried 
out within the linear response theory using the Quantum- 
ESPRESSO package [35|,[36|. Our tests show no magnetic 
moment in relevant B-rich phases allowing us to do the 
EA, phonon, and e-ph simulations without spin polariza- 
tion. 

Figure 1(a) summarizes calculated T = K forma- 
tion energies of the considered Fe^Bi^ ordered struc- 
tures with the convex hull drawn (in cyan) through the 
known oP8-FeB and tI12-Fc2B ground states; other rel- 
evant oIlO, oPIO, mP8, hP3, hP6, 0P6, oP12, 0S8, tI16, 
tI32 and cF116 structures correspond to the CTB4, FeB4 
(proposed), FCB3 (proposed), AIB2, RCB2, RuB2, FeB2 
(proposed), CrB, MoB, Ni 3 P, and Cr 2 3C6 prototypes, re- 
spectively. Phases with x > 0.5 show an expected order- 
ing, with tI12-Fe2B being stable and Fe3B and Fe23B@ be- 
ing metastable by less than 20 meV/atom. For x < 0.5, 
we find a set of phases that are below or close to the 
cv-Bo-oP8-FeB tie-line to be viable ground state candi- 
dates. We discuss their relative stability using structural 
and electronic density of states (DOS) information shown 
in Figs. 2 and 3. 

The similarity of the local coordinations in 0P8, 0S8, 



and tI16 at 1:1 composition was discussed previously in 
Ref. 37 1. In tI16 the B chains extend in two directions, a 



feature that could differentiate the structure's mechanical 
response to external load from the behaviour of the other 
two polymorphs. Fig. 1(d) reflects the difference in the 
vibrational properties of 0P8 (0S8) and tI16 and could 
explain why tI16-FeB, marginally the most stable phase 
at T = K in our calculations, has never been observed. 
The electronic DOS of FeB in the three configurations 
are rather similar (Fig. 3): they are all metallic and 
magnetic and have bonding p-d hybridized states in the 
—6 — 3 eV energy range. 

At 1:2 composition, all known metal borides stable un- 
der normal conditions are composed of 2D boron layers 
that are flat in hP3, armchair in 0P6, zigzag in hP6, 
or mixed in hR18 [19|. A detailed rigid band approxi- 
mation study of the TMB2 phases linked the distortion 
of the B layers to population of antibonding TM-TM 
and TM-B orbitals in hP3-TMB 2 with high d-electron 
count 0. The projected DOS in hP3-FeB 2 (Fig. 3) 
shows a mismatch in the maxima of the filled B and 
Fe states in the —7-— 3 eV range and a high DOS at 
the Fermi level resulting in a magnetic moment of 0.26 
/is/atom. The magnetisation energy of 11 meV/atom is 
insufficient to stabilize the hP3-FeB phase which leaves it 
200 mcV/atom above the a-BooP8-FeB tie-line. Puck- 
ering of the B layers in non-magnetic 0P6 and hP6 proves 
to be a more favorable way of reducing the high DOS at 
the Fermi level: Fig. 3 shows a higher hybridization of 
the B-p and Fe-d states with the antibonding p — d states 
now lying just above Ep. The net result of the more 
bimodal shape of the DOS is a 217 meV/atom gain in 
stability. Even more dramatic structural and electronic 
changes take place in oP12-FcB2 discovered in our EA 
search. The disintegration of the B layers opens up a 
~ 0.5 eV band gap (likely underestimated in our semilo- 
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FIG. 3: (Color online) Calculated density of states in selected 
iron borides (the lower five compounds are non- magnetic). 



cal DFT treatment [38j ) and leads to an additional 33 
me V/ atom gain in stability. This finding rules out the 
existence of hP3-FeB 2 that has been a subject of contro- 
versy [H; oP12-FcB 2 is below the a-B^oP8-FeB tic-line 
by over 30 meV/atom (Fig. 1(c)) in the whole T range 
and should be synthesizablc. 

At 1:3 composition, the EA suggests a new mP8 
phase [13 that breaks the a-BooP8-FeB tie-line by 15 
meV/atom; however it is found to be metastablc w.r.t. 
a-B and oP12-FcB2 at all temperatures. 

At 1:4 composition, the observed proximity of the oIlO- 
FeB4 phase to the a-BooP8-FcB tie-line is intriguing as 
there have been explicit references to unsuccessful at- 
tempts to synthesize this phase [3). Unexpectedly, a 
phonon dispersion calculation showed dynamical insta- 
bility of oI10-FeB4, with imaginary frequencies reach- 
ing 208z cm" 1 for a T-point phonon in the conventional 
10-atom unit cell. Using the phonon eigenvector that 
skews the rectangle B building units in the x-y plane 
into parallelograms we have constructed a new structure 
type, 0PIO-FSB4. The considerable energy gain of 28 
meV/atom and no imaginary frequencies in the phonon 
spectrum make 0PIO-F8B4 thcrmodynamically and dy- 
namically stable in the considered temperature range rel- 
ative to known phases [Fig. 1(b)]; with phonon correc- 
tions included 0PIO-FCB4 lies 3 meV/atom above the a- 
BooP12-FeB 2 tie-line at T = K but 10 meV/atom be- 
low the tie-line at T = 900 K. Additional no-symmetry 
relaxations of distorted oIlO and oPIO supercells with 
10, 20, and 40 atoms have consistently produced oPlO- 
FeB4 as the most stable configuration. The EA search 
has also shown that Fe 2 Bs cells converge to oP10-FeB4 
while FC3B12 cells evolve into a new mS30-FcB4 phase 
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FIG. 4: (Color online) Top: total and projected phonon den- 
sity of states (PHDOS) in oP10-FeB 4 . Bottom: Eliashberg 
function and the strength of the electron-phonon coupling. 



metastable by 6 meV/atom. The tests seemed necessary 
due to a counterintuitive evolution of the DOS in the oIlO 
to oPIO transformation: the Fermi level in 0PIO-F0B4 
catches the edge of the antibonding Px,y-dx 2 -y 2 peak re- 
sulting in a high n(Ep) = 1.0 states/(eV spin f.u.); the 
feature is unusual as stable compounds tend to have the 
Fermi level lying in the pseudogap . 

The naturally electron-doped 0PIO-FCB4 candidate 
material with strong covalent bonds is next analyzed for 
superconducting features. We use the linear response 
theory and fine k- and g-meshes 27, 36 



to calculate the 
phonon DOS (PHDOS), Eliashberg function (a 2 F(uj)), 
and strength of the e-ph coupling (A(w)). The phonon 
spectrum in Fig. 4 can be divided into three regions 
with mixed Fe-B modes (0-320 cm -1 ), B modes with a 
relatively flat PHDOS (320-740 cm" 1 ), and B modes in- 
volving in-plane optical distortions of B parallelograms 
(860-920 cm" 1 ). The Eliashberg function integrates to 
a large Atot = 0.80 and gives the logarithmic average 
= 430 cm" 1 . While key contributions to A to t ~ 0.8 
in CaC6 and MgB 2 come from the low-frequency Ca 
modes (co < 150 cm" 1 ) [4l[ and the high-frequency B 
modes (500 cm" 1 < uj < 560 cm" 1 ) [42|, respectively, 
nearly 60% of Atot m oP10-FeB4 is generated by the 
mixed Fe-B modes in the 160-300 cm" 1 range, (w), 
in oP10-FeB4 is found to be much closer to the MgB 2 
value of ~ 450 cm -1 [4^], rather than the CaCg value 
of — 200 cm -1 [Il|. Using the Allen-Dynes formula 
and typical \x* of 0.14-0.10 we estimate the T c in oPlO- 
F3B4 to be 15-20 K. The compound has two 3D Fermi 
surfaces centered at T and R points and the T c may be 
further enhanced by the multiband effect. Because of 
the large gradient of the DOS near the Fermi level the 
superconducting properties may be strongly affected by 
the presence of vacancies or impurities. Although oPlO- 
FeB4 is found to have neither ferro- nor antiferromagnetic 
moment, spin fluctuations 23j could play a critical role in 
the pairing mechanism and should be examined carefully 
using input from experiment. 
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In summary, our search for new compounds in the 
common Fe-B system demonstrates the necessity to go 
beyond standard structure types: The EA-driven uncon- 
strained structural optimization has uncovered a set of 
viable ground states, with new oP10-FeB4 and oP12- 
FeB2 shown to be thermodynamically stable by over 
25 meV/atom relative to the known a-B and oP8-FeB. 
To the best of our knowledge, the identified boron-rich 
phases have been never observed before and their dis- 
covery may require finding suitable kinetic routes. The 
presented analysis of the structural and electronic proper- 
ties shows how the phases stabilize and what new physics 
they are expected to exhibit if synthesized, (i) 0PIO-FCB4 
could become yet another exception to Matthias' rules 
[44} that recommend staying away from magnetic ele- 
ments when designing new superconductors. This com- 
pound has a high DOS at the Fermi level leading to 
A = 0.80 and a surprisingly high T c of 15-20 K. (ii) oP12- 
FeB2 is predicted to be the first metal diboride semicon- 
ductor with a ~ 0.5 eV band gap (time-dependent DFT 
or Green's functions techniques will likely give a larger 
value 38]). (iii) The proposed materials may also exhibit 
appealing mechanical properties as hardness tends to be 
higher for M x Bi_ x with x < 0.5 [22TJ . 
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Library of structure types We have considered over 40 
commonly seen M-B and M-C structure types listed 
in the ICSD. The list below is ordered by composi- 
tion and specifies the Pearson symbol and the proto- 
type for each entry. 1:12: cF52-UBi 2 ; 1:7: oI64-MgB 7 ; 
3:20: oS46-Na 3 B 2 o; 1:6: cP7-CaB 6 ; 3:14: tI160-Li 3 Bi 4 ; 
1:4: tP20-UB 4 , oI10-CrB 4 , mS10-MnB 4 ; 1:3: tP20-LiB 3 , 
hP16-Moo. 8 B 3 ; 2:5: hR21-Mo 2 B 4 . 65 ; 1:2: hP3-AlB 2 , 
oP6-RuB 2 , hR18-MoB 2 , oP12-PbCl 2 , hP6-ReB 2 , hP6- 
MoS 2 ; 2:3: oS20-V 2 B 3 , hP10-Ru 2 B 3 ; 3:4: oI14-Ta 3 B 4 ; 
5:6: oS22-V 5 B 6 ; 1:1: oP8-FeB-b, oS8-PtB . 6 7, tI16- 
MoB, 0S8-TII, oP8-TiSi, hP2-WC; 5:4: hP18-Rh 5 B 4 ; 
11:8: oP38-RunB 8 ; 3:2: tP10-U 3 Si 2 , aP30-Rh 9 B 5 . 5 ; 5:3: 
tI32-Cr 5 B 3 ; 2:1: tI12-CuAl 2 ; 7:3: hP20-Th 7 Fe 3 , oP40- 
Mn 7 C 3 ; 5:2: mS28-Mn 5 C 2 ; 3:1: oP16-Fe 3 C, oS16-Re 3 B, 
tI32-Ni 3 P; 23:6: cF116-Cr 23 C 6 ; 4:1: tP10-Be 4 B. TableU 
summarizes the lattice parameters and Wyckoff positions 
of the fully relaxed unit cells of the relevant ICSD and 
proposed phases. 

Evolutionary algorithm (EA ) search Module for Ab Ini- 
tio Structure Evolution (MAISE) [l[ takes advantage of 
the best EA strategies reported in the recent years @- 
5] as well as our extensive tests to achieve fastest con- 
vergence for metal borides. The module creates input 
for VASP and gets back fully relaxed structures with 
the corresponding ab initio enthalpies. Starting sets 
of N — 12 structures are generated randomly. The 
population is then evolved through crossover (75%) and 
crossover plus mutation (25%); the mutation consists of 
shifts and swaps of the atoms as well as distortions of 
the unit cell. The best N members are selected from 
the pool of N old and N new members; the probabil- 
ity of survival is proportional to the fitness depending 
on the structure's enthalpy. Special care is taken to de- 
tect and eliminate similar structures in order to avoid 
stagnation. Fig. l(a,b) demonstrates that diversity is 
maintained throughout the search (grey points) while the 
fittest member is always kept in the population (green 
points). Fig. 1(c) shows the lowest relative formation 
enthalpies obtained in three runs with different random 
seeds for each composition. Each run consisted of 50-100 
generations, i.e., 600-1200 structural relaxations. Typi- 
cally, all three runs converged to the same structure in 
no more than 60 generations. 

Accuracy of the DFT approximations The 500 eV en- 
ergy cutoff and dense fc-meshes ensure numerical con- 
vergence of the formation energy differences to within 
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FIG. 1: a,b) Typical enthalpy distributions for competing 
structures during an EA-driven ground state search. In each 
generation 12 new members are fully relaxed with VASP; the 
best member shown in green is kept in the population. For 
1:4 and 1:2 Fe-B compositions the evolutionary search finds 
new structure types, oPIO and oP12, significantly lower in 
enthalpy than the known ones from the ICSD. c) Relative for- 
mation enthalpies (referenced to the a-B<->-oP8-FeB tie-line) 
of the most stable phases for each unit cell size found in the 
EA search; compounds with negative AHf orm are thermody- 
namically stable w.r.t. the known a-B and oP8-FeB. 



a few meV/atom but the systematic errors due to the 
approximated form of the xc functionals for the com- 
plex Fe-B system can be expected to be larger. Fig- 
ure 1 shows dependence of the relative stabilities for 
predicted phases on the choice of commonly used DFT 
flavours: the Perdew-Zunger parameterization within the 
local density approximation (LDA) Q and the Perdew- 
Burke-Ernzerhof (PBE) 0, Perdew-Wang (PW91) 0], 
and PBEsol [9] parameterizations within the generalized 
gradient approximation (GGA). 

We observe that the relative stability of the boron-rich 
phases varies by over 30 meV/atom when simulated in 
GGA-PBE and LDA but in all the treatments the pre- 
dicted FeB 2 and FeB 4 compounds are stable w.r.t. the 
a-B-H>oP8-FeB tie-line. The systematic errors do not al- 
low us to conclusively resolve the close in energy 0P6- 
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FIG. 2: E(V) curves for the competing crystal structures 
oP8, oS8, and tI16 at 1:1 composition relative to oP8 (upper 
panel) and the predicted structures at 1:2 and 1:4 composition 
(lower panel) relative to the a-B-ooP8-FeB tic-line. 



FeB 2 and hP6-FeB 2 or tI16-FeB, oP8-FeB and oS8-FeB. 
However, oP10-FeB4 is consistently below oI10-FeB4 by 
~ 25 meV/atom and hP3-FeB2 is consistently above the 
a-BooP8-FeB tie-line by ~ 200 meV/atom. The ther- 
modynamic corrections are found to be not sensitive to 
the choice of the DFT approximation and show similar 
T-dependences in the LDA and GGA-PBE. 

Convergence of the e-ph calculations In order to check 
the convergence of A we perform the calculations for two 
sets of grids: (i) the 2 x 2 x 4 q- and 6x6x12 fc-mcshcs 
for the dynamical matrix and the 12 x 12 x 24 k— mesh 
for the e-ph coupling (ii) the 3 x 3 x 6 q- and 9 x 9 x 18 
fc-meshes for the dynamical matrix and the 18 x 18 x 36 
fc-mesh for the e-ph coupling. For the 0.02 Ry smearing 
in the /c-point integration we find nearly identical values 
of A = 0.80 (note that in these calculations the smearing 
leads to a lower value of n(Ep) « 0.8 states/(eV spin 
f.u.) due to a large gradient of the DOS near the Fermi 
level). We also find negligible changes (below 0.4 cm" 1 ) 
in T-point frequencies when the fc-mesh is increased from 
9 x 9 x 18 to 12 x 12 x 24 and 15 x 15 x 30. 
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TABLE I: Space groups, lattice parameters, and Wyckoff-positions of fully relaxed Fe-B structures. 
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